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Abstract 

In this paper, we investigate the spectral instability of periodic traveling wave so- 
lutions of the generalized Korteweg-de Vries equation to long wavelength transverse 
perturbations in the generalized Kadomtsev-Petviashvili equation. By analyzing high 
and low frequency limits of the appropriate periodic Evans function, we derive an ori- 
entation index which yields sufficient conditions for such an instability to occur. This 
index is geometric in nature and applies to arbitrary periodic traveling waves with minor 
smoothness and convexity assumptions on the nonlinearity. Using the integrable struc- 
ture of the ordinary differential equation governing the traveling wave profiles, we are 
then able to calculate the resulting orientation index for the elliptic function solutions 
of the Korteweg-de Vries and modified Korteweg-de Vries equations. 

1 Introduction 

The Korteweg-de Vries equation 

(1.1) Ut = Uxxx + UU^ 

often arises as an model for one-dimensional long wavelength surface waves propagating in 
weakly nonlinear dispersive media, as well as the evolution of weakly nonlinear ion acoustic 
waves in plasmas |TRRj . When the assumption that the wave is purely one-dimensional 
is relaxed to allow for weak dependence in a transverse direction one is led to a variety 
of multidimensional generalizations of the KdV. One of the most well studied weakly two- 
dimensional variations of the KdV is the Kadomtsev-Petviashvili (KP) equation [KP] given 
by 

(1.2) {Ut - UxXX - UUx)^ + (TUyy = 0, 
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where the constant a = ±1 differentiates between equations with positive {a = +1) and 
negative (<t = —1) dispersion and the choice of a depends on the exact physical phenomenon 
being describecOl- For instance, if a = +1 (ll.2p is referred to as the KP-I equation, which 
can be used to model waves in thin films with high surface tension, while (jl.2p is called 
the KP-II equation in when a = —1, which can be used to model water waves with small 
surface tension. Various other physical applications utilize equations of the form (jl.2p . such 
as the modeling of small amplitude internal waves and in the study of unmagnetized dusty 
plasmas with variable dust charge |PJ] . 

In many applications, appropriate scaling in the physical parameters introduces a pa- 
rameter a > in the nonlinearity yielding a governing equation of the form 

Ut = Uxxx + CtUUx- 

In neighborhoods of parameter space where a = one is forced to consider higher order 
expansions in the nonlinearity, the most natural being of the form 

Ut = Uxxx + Pu^Ux 

where /3 7^ 0. The choice of the sign of (3 must be made depending on the particular physical 
situation being studied. In particular, the choice of j3 clearly determines the structure of 
the stationary homoclinic/heteroclinic solutions and hence the case of /3 > and /3 < 
define quite distinct dynamics. In this case the modified KP equation 

{ut - Uxxx - Pu^Ux)^ + auyy = 

arises naturally as a weakly two-dimensional generalization. In the context of dusty plasmas 
with variable dust charge, the mKP equation can be derived near the "critical" density case 
(see [P Jj ) for details. 

As our theory will not depend on the explicit form of the nonlinearity, and to encompass 
as many physical applications as possible, we will most often work with the generalized KdV 
(gKdV) equation 

(1-3) Ut = Uxxx + fiu)x 

and its weakly two-dimensional variation the generalized KP equation 

(1-4) (ut - Uxxx - f{u)x)x + ^Uyy = 

where the nonlinearity / is sufficiently smooth and satisfies general convexity assumptions. 
For such nonlinearities, the gKdV equation admits asymptotically constant traveling solu- 
tions, known as solitary waves, as well as traveling waves which are spatially periodic. It 
is the latter case we consider here. As a solution u{x, t) of (II. 3p is clearly a y-independent 

^Notice that one can always rescale the y variable to force a to be any (non-zero) real number. However, 
for convenience, we will always assume that a — ±1. 
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solution of ()1.4p . it seems natural to question the stability of such a solution to perturba- 
tions which have a nontrivial dependence on the transverse (y) direction. Such a transverse 
instability analysis is the subject of the current paper: in particular, we study the spectral 
stability of a stable y-independent spatially periodic traveling wave solution of the gKdV 
equation to perturbations which are co-periodicH in x with low frequency oscillations in the 
transverse direction. To this end, we will loosely follow the general Evans function approach 
of Gardner [Glj with additional aspects of analysis from the more recent work of Johnson 

The transverse instability of solitary waves of the KdV in the KP equation was first 
conducted by Kadomtsev and Petviashvili [KP] , where it was found that such solutions are 
stable to transverse perturbations in the case of negative dispersion, while they are unstable 
to long wavelength transverse perturbations in the case of positive dispersion (even though 
they are stable in the corresponding one-dimensional problem). Moreover, in |APS] it was 
shown that a short wavelength cutoff for instability exists for the positive dispersion case 
and the dominate mode of instability was identified. Other authors have devised various 
techniques to demonstrate the transverse instability of KdV solitary waves in the KP-I 
equation: for example, see the pioneering work of Zakharov |Za] . where the author utilizes 
the integrability of the KP-I equation via the inverse scattering transform, and the recent 
work of Rousset and Tzvetkov |RT2| , where the authors use general PDE technique^. The 
transverse instability of a one-dimensionally stable gKdV solitary wave in the corresponding 
gKP equation has recently been considered in |KTN] using perturbation analysis similar 
to that of Kadomtsev and Petviashvili. In particular, multiple scale analysis was used to 
derive an evolution equation for the wave velocity to describe the slow-time response of 
the solitary wave in response to the long wavelength transverse perturbations. The authors 
conclude that for positive dispersion the solitary waves of the gKdV equation are always 
unstable to long wavelength transverse perturbations in the gKP equation. Moreover, it 
was found that for some nonlinearities the solitary waves may in fact be unstable in the 
case of negative dispersion. 

In the case where the background solution of the gKdV is spatially periodic, there do 
not seem to be any results concerning the transverse instability in the gKP equation and, 
moreover, the stability of such solutions in general is much less understood. This refiects 
the fact that the spectrum of the corresponding linearized operators is purely continuous, 
and hence it seems more difficult for nonlinear periodic waves of the gKdV to be stable than 
their solitary wave counterparts. Moreover, the periodic waves of the gKdV in general have 
a much more rich structure than the solitary waves: even in the case of power-law nonlin- 
earities it is not possible to write down a general elementary representative for all periodic 
solutions of the gKdV, which stands in contrast to the solitary wave theory. Nevertheless, 
there has been much study recently into the one-dimensional stability of such solutions 

^Notice by Floquet theory, one-dimensional spectral instability to co-periodic perturbations implies spec- 
tral instability to localized perturbations. 

^The methods of Rousset and Tzvetkov provide nonlinear transverse instability of the KdV solitary wave, 
and has the advantage over that of Zakharov of generalizing to the full water wave problem in the presence 
of surface tension. 
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under the gKdV flow with results ranging from spectral stability to localized perturbations 
(see [BDj . |Br J] . and |HK| ) to nonlinear (orbital) stability to co-periodic perturbations (see 
[Br JK] ■ |DKj . and [Jl|). We also want to stress the fact that the periodic solutions are 
also physically relevant as they are used to model nonlinear wave trains: such patterns are 
prevalent in a variety of applications and their instabilities have been studied extensively 
in the literature (see for example the classic works of Benjamin [Be2], Benjamin and Feir 
[BF] . Lighthill [L], and Whitham[W].). 

The goal of this paper is to perform a transverse instability study in the case of spatially 
periodic traveling wave solutions which are stable to perturbations in the x direction of uni- 
directional propagation. To this end, we will develop a somewhat nonstandard orientation 
index which detects instability of a T-periodic traveling wave of the gKdV to perturbations 
in the gKP equation which are T-periodic in the x-direction with low frequency oscillations 
in the transverse spatial variable y. Throughout this paper, we will refer to such pertur- 
bations as long wavelength transverse perturbations. This is accomplished by studying 
the behavior of the corresponding periodic Evans function in both the high (real spectral) 
frequency and low (transverse) frequency limits when the wave number of the transverse 
perturbation is small but non-zero. As we will see, the high frequency analysis is somewhat 
delicate due to a degeneracy in the KP equation. In particular, it is seen that one must take 
into account not only the higher order effects of the dependence of the limiting asymptotic 
ordinary differential equation on the background solution as in |MaZ3| and [PZ] . but also 
the inherent averaging/cancellation effects due to the periodicity. Such a result seems new 
to the literature, and relies on the use of block-triangularizing transformations as used in 
[HLZj . Moreover, in the appendix we give a slight simplification of the of the tracking lem- 
mas utilized in [MaZ3] . |PZ] . and [HLZj by formulating the result in terms of conjugating 
transformations rather than invariant graphs. This view point seems more closely related 
to the style and needs of current research (see for example the related treatments in |Z6] 
and [NZ]). As a result, we will find that the limiting behavior of the sign of the periodic 
Evans function for large (real) spectral frequency is governed precisely by the sign of the 
dispersion parameter a. Thus, the stability of our periodic traveling waves depends directly 
on the case of dispersion chosen in the gKP equation, and hence on the exact physical 
phenomenon being described. 

The low frequency analysis utilizes matrix perturbation theory and the methods of |BrJ] 
and [J2]. In particular, we will explicitly construct four linearly independent stationary 
solutions of the linearization of the governing PDE: three of these will be given to us by 
variations in the traveling wave parameters, which we will use to parameterize the traveling 
wave solutions of ()1.3p . while the fourth can be constructed using standard techniques. 
Using variation of parameters then, we can compute the leading order variation of these four 
functions in the transverse wave number k for |A;| <^ 1, and hence can determine the leading 
order variation of the periodic Evans function in the transverse wave number at the origin 
in the spectral plane. This leading order variation is expressed as a Jacobian determinant 
relating to the ability to parameterize near by periodic waves of fixed wavespeeds by the 
conserved quantities of the gKdV flow. It follows then that if this determinant has the 
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opposite sign of that of the high frequency hmit (the parameter a), we immediately have a 
spectral instability to long wavelength transverse perturbations in the gKP equation. Notice 
this approach is somewhat different than that of [J2], where transverse instability in the 
generalized Zakharov-Kuznetsov equation was established by finding sufficient conditions 
for eigenvalues bifurcating from the origin (in k) to enter the unstable half plane. In that 
case, the low frequency analysis involved considering variations of elements of the kernel of 
the linearized operator in both the spectral variable and in the transverse wave number. As 
a result, the results in [J 2] concern spectral instability to long wavelength perturbations in 
both the direction of propagation of the background solution and the transverse direction. 
The orientation index derived in our case thus seems to be nonstandard, in the sense that 
it detects instabilities to perturbations which have bounded period in the x-direction and 
admit slow modulations in the transverse direction. 

The outline of this paper is as follows. In section 2, we review the basic properties 
of the periodic traveling wave solutions of the gKdV equation ()1.3p . In particular, we 
will discuss a parametrization of the traveling wave solutions of (jl.Sp which will be useful 
throughout our analysis. In section 3, we conduct our transverse instability analysis by 
first conducting a high frequency analysis of the associated periodic Evans function and 
then conducting the corresponding low frequency analysis. As a result, we will have an 
instability index which guarantees that a periodic traveling wave solution of the gKdV is 
spectrally unstable to long wavelength transverse perturbations in the gKP equation. This 
index can be calculated exactly in several model cases and the relevant results will be given. 
In section 4, we end with some closing remarks and in the appendix we present a proof of 
the tracking lemma utilized in the high frequency analysis and also outline the proof of a 
formula crucial to the low frequency analysis. 

2 Properties of Periodic Traveling GKdV Waves 

In this section, we review some of the basic properties of the periodic traveling wave solutions 
of the gKdV equation. For more details, see [BrJ] or [Jlj . For each c > 0, a traveling wave 
with speed c is a solution of the ordinary differential equation 

(2.1) Uxxx + fiu)x - CUx = 0, 

i.e. they are stationary solutions of (jl.Sp in the moving coordinate frame defined hy x + ct. 
This equation is clearly Hamiltonian and hence we can reduce it to quadrature. Indeed, by 
integrating (|2.ip twice we see that a traveling wave profile of the gKdV must satisfy the 
nonlinear oscillator equation 

2 

(2.2) !^ = E + au + ^u^ -F{u), 

where F is an antiderivative of the nonlinearity / satisfying F{0) = and a and E are 
constants of integration. Thus, the traveling waves form a four parameter family of solutions 
of (jl.Sp described by the constants a, E, and c together with a fourth constant of integration 
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corresponding to a translation mode: this translation direction is simply inherited from the 
translation invariance of (jl.3p and hence can be modded out. It follows that on open subsets 
of M.^ = (a, E, c) equation (I2.2p admits a periodic orbit. Moreover, the boundary of these 
open subsets correspond to solutions which decay asymptotically at infinity and hence can 
be identified with the solitary wave solutions. In particular, notice that in order for ()2.2p to 
admit a solitary wave solution, the constant E must be fixed by the prescribed boundary 
conditions and we must have a = 0. It follows the solitary waves form a codimension two 
subset of the family of traveling waves. 

In general, equation ()1.3p admits three conserved quantities. In order to define these, 

let 

V{u; a, c) = F{u) — au — -v? 

be the effective potential arising in the nonlinear oscillator equation (j2.2p . Throughout 
this paper, we will assume the roots u± of the equation E = V{u; a, c) are simple, satisfy 
< and that V{u;a,c) < E for u E («_,«+). As a consequence, u± are functions 
of the traveling wave parameters a, E, and c and, without loss of generality, we can set 
u{0) = n_. It follows that we can express the period of a periodic solution of ()2.ip via the 
formula 

T = T{a,E,c) = V2 ^ 

Ju. E — V{u;a,c) 

By a standard procedure, the above integral can be regularized at the square root branch 
points and hence represents a function of a, E, and c. Similarly, the conserved quantities 
of the gKdV flow can be represented as 



M{a,E,c) = [ u{x) dx = 2 [ 

Jo Ju 

P{a,E,c)= [ u^{x)dx = 2 [ 

Jo Ju 



u du 



^2{E-V{u-a,c)) 
^ V? du 



yj2{E -V{u;a,c)) 

i/(a, E,c)= — - F(u) ]dx = 2 / — , \ ^ ^ J — \_ L 

^ Wo V 2 ^7 Ju. y/2{E-V{u;a,c)) 

representing the mass, momentum, and Hamiltonian, respectively. As above, these integrals 
can be regularized at the branch points and hence represent functions of the traveling 
wave parameters. As we will see, the gradients of the period and mass of the solution u will 
play a very large role in this paper. However, as pointed out in |Br J) . when E ^ gradients 
in the period can be interchanged for gradients of the conserved quantities via the relation 



EVa,E,cT + aVa,E,cM + ^Va,E,cP + Va,E,cH 



■ 0. 



where Va,E,c = {da, Qe, dc) ■ Thus, all gradients involved in the results of this paper can 
be expressed completely in terms of the gradients of the conserved quantities of the gKdV 
flow, which seems to be desired from a physical point of view. 



3 TRANSVERSE INSTABILITY ANALYSIS 



7 



We now discuss our parametrization of the family of periodic traveling wave solutions 
of (jl.Sp more carefully. A major technical assumption throughout this paper is that the 
period and mass provide good local coordinates for the periodic traveling waves of fixed 
wave speed c > 0. More precisely, given a periodic traveling wave «(•; oq, -Eq, cq) of (jl.Sp 
with Co > we assume the map 



have a unique inverse in a neighborhood of (aoj-Eo) £ K^, which is clearly equivalent 
with the non- vanishing of the Jacobian determinant 



at the point (ao, Eq, cq). As we will see, the sign of this Jacobian controls the low frequency 
analysis presented in this paper. It is worth mentioning that while this may seem like a 
rather obscure requirement, this Jacobian has already shown to be important in the stability 
theory of periodic traveling waves of the gKdV: see for example |Br J] . |Br JK] . |J1) . and |J2) . 
In particular, this Jacobian been computed in [BrJKj for several power-law nonlinearities 
and, in these cases, has been shown to be generically non-zero. Moreover, such a non- 
degeneracy condition should not be surprising: a similar non-degeneracy condition must 
often be enforced in the stability theory for solitary waves (see |Boj . |Be2] . and |P W] ) . 

3 Transverse Instability Analysis 

We now begin our stability analysis. Let u = u{-;a, E,c) be a T = T(a, c)-periodic 
traveling wave solution of (jl.3p . Moreover, we assume that u is a stable solution of the 
one-dimensional gKdV equatior|f|. As noted in the introduction, it is clear then that u is 
a y-independent solution of the generalized KP equation (jl.4p for either cr = ±1. We are 
interested in the spectral stability of u as a solution of (II. 4p to small perturbations. To this 
end, consider a small perturbation of u of the form 

i{j{x, y, t) = u{x) + £v{x, y, t) + ©(e^), |e| < 1 

where v{-,y,t) £ L^{R) for each {y,t) £ and v{x,-,t) G L°°(M) for each {x,t) G R"^. 
Forcing ip to solve the traveling gKP equation 

(3.1) {Ut - Uxxx - f{u)x + CUx)^ + (TUyy = 0, 

yields a hierarchy of consistency conditions. The 0{e^) equation clearly holds since u solves 
(13. ip . and the 0{e^) equation reads as 



{a,E) ^ {T{a,E,co),M{a,E,co)) 




dx {dt + dxC[u]) V + avyy = 



*Else the issue of transverse instability is of no interest, as instabilities to unidirectional perturbations 
will prevent stability to higher dimensional perturbations. 
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where C[u\ = —d^ — f'{u) + c is a periodic Hill operator. As this linearized equation is 
autonomous in both time and the spatial variable y, we may seek separated solutions of the 
form 

v{x,y,t) = e->'^-^'''yv{x) 
where G C, /c G M, and v £ L^(IR). This leads one to the (generalized) spectral problem 

(3.2) {dlC[u] - ak^) v = fid^v 

considered on the real Hilbert space L^(IR). We refer to the background solution u as being 
spectrally stable in L'^{R) if (fO) has no L'^{R) spectrum witlil 5R(/i) / for any k£R. 

Since the coefficients of the differential operator C[u] are T-periodic, as they depend on 
the background solution u, standard results in Floquet theory implies the L^ spectrum of 
()3.2p is purely continuous, and consists entirely of L°°(M) eigenvalues. Indeed, the fact that 
(j3.2p can have no L^(]R) eigenvalues is clear: writing (13. 2p a first order system of the form 

and letting ^>(x; n, k) be a matrix solution satisfying the initial condition <I?(0; n,k) =1 for 
all (/z, fc) G C X R, we define the monodromy operator, or the period map, to be 

M(^) :=<I>(T;^,A;). 

Notice that given any vector solution Y of ()3.2p . the monodromy operator is a matrix such 
that 

M(^)Y(x; n, k) = Y{x + T; /i, k) 

for all {x, fj,,k) G M x C x R. Assuming now for simplicity that Y is an eigenvector of M(^) 
with eigenvalue A, we clearly have that 

Y{NT; fi, k) = M(^)^Y(0; ^, k) = A^Y(0, ^, A;). 

Thus, if Y{x;fi,k) decays as x — )• oo it must become unbounded as x — >• — oo. Thus the 
best we can hope for is for Y(x; k) to remain bounded on M, which corresponds in this 
example to A G 5^, i.e. |A| = 1. For more details, see [Hj for example. 

Following Gardner (see |G1] and |G2] ). we define the periodic Evans function for our 
problem to be 

D (/i, /c. A) = det (M(/i, k) - AI) , (/i, /c. A) G C x M x C. 

The complex constant A is called the Floquet multiplier and is related to the class of 
admissible perturbations in (j3.2p . In particular, notice that A = 1 corresponds to T-periodic 
perturbations of the background solution u. Clearly, Z)(/U, k, A) is an entire function of ^ and 

^Usually, one defines spectral stability as the absence of spectrum with positive real part. However, in 
our case the spectrum is symmetric about the imaginary axis and hence spectral stability is equivalent with 
the spectrum being confined to the imaginary axis. 
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k for each fixed A € C since the coefficient matrix H(x,/i, /;;) depends as such on /i and k. 
This allows an analytical characterization of the L^(M) spectrum of (j3.2p : the generalized 
spectral problem (13.20 has a non-trivial bounded solution for a given A; G M if and only if 
there exists a k G M such that 

D{^i,k,e''') = 0. 

In particular, D{fi,k,l) = if and only if ()3.2p has a non-trivial bounded solution in 
Lpg;,([0, r]) for a given /c G M. Moreover, the following property will be useful in the low 
frequency analysis conducted later in the paper. 

Lemma 3.1. The function D{fj,, k, A) is an even function of both /i and k. 

Proof. Since the spectral problem (c?^/^[?x] + fj/c^) v = fidxV is invariant under the transfor- 
mation k I—)- —k, it follows that D{fj,,k,\) is an even function of k. To analyze the parity 
in n, we write 

Dii^i, k, A) = det (M(^, k) - AI) 

= A'' + a{fi, k)X^ + k)X^ + c{fi, k)\ + 1 

where a{ii,k) = -tr(M(^,A;)) and b{fj.,k) = ^ ( tr(M(/i, /c)^) - tr(M(/i, /c))^) . Since the 
spectral problem is invariant under the transformation x — t- —x and fi — t- —fi, it follows that 
the matricies M(/u, A;) and M(— /u,A;)~^ are similar and hence a direct calculation yields 

D{fi, k, A) = A^ det ( M(-/i, k) - \l 



A 

= A^ + c(-^, A;)A^ + b{-fi, k)\^ + a(-^, A;)A + 1. 
It follows that c(//, A;) = a{—fi, k) and b{n, k) = 6(— /x, k), and hence 

D{^l, k,l) = 2 + a(^, k) + a(-M, k) + Y^"^'^) . 

Thus, /c, 1) is an even function of /x. □ 

The goal of our analysis is to provide sufficient conditions to ensure that when < \k\ <^ 
1 the function /j, i— >■ -D(/U, k, 1) has a non-zero real root, corresponding to an exponential 
instability of the underlying wave. To this end, we derive an orientation index by comparing 
the high frequency and low frequency (in /x) asymptotics of the function D(fj,,k,l). In 
particular, we will see that the sign of k, 1) for large real /i equals the sign of a. Thus, 
if the quantity -D(0, k, 1) has the opposite sign of a, we can infer the existence of a /x* G M"*" 
such that D{fj,*,k, 1) = which implies exponential instability of the background solution. 
We begin by analyzing the high frequency behavior of the periodic Evans function. 
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3.1 High Frequency Limit 



In this section, we study the large real fi behavior of the periodic Evans function D(fi, k, 1) 
when k ^ 0. To begin, rescale (j3.2p with the change of variables x = \fi\^^^x to obtain (after 
dropping the tilde's) the spectral problem 



(3.3) {-dt - \fi\-'/'dl (/'(u) + c) - ae\^^\-'/^ 

This can be rewritten as a first order system of the form 



a 



V. 



(3.4) 



/ 1 \ 

10 

1 

V -1 y 



w + 



/ 











w. 



where Ai = — 2/"(u)n^, A2 = —f'{u) + c and 



X 



'2/3 



b(m) 



-4/3 



On a heuristic level then, we expect that the monodromy operator for ^ S> 1 will behave 
like 

M(^) ^ e^^^\^\'"^ 

and hence 

L»(//, fc, 1) « det (e^ol'^l'^'^ - l) 

However, the matrix Hq clearly has an eigenvalue of and hence this heuristic argument 
leads us to expect that D{^, /c, 1) — as /i — +00. From the point of view of an orientation 
index, this does not provide us with sufficient information: we must know what the limiting 
sign of the Evans function is. Thus, although we expect D{^, k, 1) vanishes in the limit as 
^ —7- 00, we must analyze the situation more closely to determine if the sign of -D(/U, fe, 1) 
has a limiting value. This is the content of the following lemma. 

Lemma 3.2. For k ^ 0, we have the high frequency limit 

lim sgn {D{fi, k, 1)) = sgn{a). 

Remark 3.3. Notice this result is somewhat unexpected due to the form of the rescaled 
equation ()3.3p since the term ak"^ enters at only second order in |^|~^/'^. Thus, this is a 
very small term compared with the 0{\^\~'^^^) terms involved. However, the following proof 
will show that, upon averaging, the periodicity of the underlying solution implies cancelation 
of the lower order effects and hence the asymptotic behavior for large ji must he determined 
by the higher order terms. 
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Proof. First, notice that by Lemma l3.1l it is enough to consider the hmit as /i — )• +00 only. 
Consider the rescaled first order system p.4p and define e = |/u|~^/^ ^ 1. We refer to the 
constant matrix Hq as the principle part and regard the matrix B as a matrix of error 
terms. We begin by diagonalizing the principle part. To this end, let A = ^ (l + i\/3) and 
notice if we define the matrix 



then a straight forward calculation yields 



/-I 


-1 


-1 


1\ 


1 


-A 


-A* 





-1 


A* 


A 





V 1 


1 


1 


0/ 



/-I 0\ 



Q^HoQ 



A 

A* 



V 0/ 



Note that 3f?A = KA* is positive, hence the real parts of the diagonal entries of the stable, 
neutral, and unstable diagonal blocks of Q~^HoQ each have a spectral gap, one from the 
other. Using the change coordinates W = QY in (13. 4p . we can consider the first order 
system 



By a brief but tedious calculation, we see that the matrix B take the block form 
/ 





/-I 















A 
















A* 







Vo 








0/ 




B 



1 \ 

3X 



0{e) 



X 



\{Ai-A2)e-x\{-XAi + X*A2)e-x\{->^*Ai + XA2)e-x \ X J 

where the upper left hand block is 3 x 3 and the lower right hand block is 1 x 1. Now, define 
the T-periodic matrix valued function 











U + e 



-Ai + A2 + ]^Ai.^ \ -Ai + y ^2 + 



-Ax + -^^2 + TT^^i.a^ i 
A* 2A* ; 
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where I4 is the standard 4x4 identity matrix and again the upper left hand block is 3 x 3. 
Then another straightforward computation implies that 



/ 0{e 



S-^BS 



Oie) 



I 



where again the upper left hand block is 3 x 3. 

Now, noticing that S' = ©(e^/^) and that S is a 0{£) perturbation of the identity, it 
follows that making the variable coefficient change of variables U = SY yields a first order 
system of the form 



Q 'HoQ + 



0{e) 



1 



U. 



In particular, we see that the resulting coefficient matrix is approximately block upper 
triangular with error of order 0(e^/^) = Od/il"^). Applying Lemma lA. II from Appendix A 
(notice that here = 1, iV = ©(e), and b = and Remark IA.3I we find that there is 

a T-periodic change of coordinates X = ZU of the form 



/3 
$ 1 



where $ = 0(e^/^) is of dimension 1x3, taking the system to an exact upper block triangular 
form with diagonal blocks 



-l + 0{e) 



A 

A* 



+ 0(e), and ^Ai,^e + Q^iAi,, - akA + 0{t'l^) 



Finally, by the block-triangular form plus periodicity of the coordinate changes we may 
compute the periodic Evans function as the product of the periodic Evans function of the 
diagonal blocks of this transformed system, integrated over a period T := r|//|^/^ going to 
infinity. The contribution from the stable block i^ approximately e"'^'^'^^^ — 1, so has sign 
— 1. Similarly, the unstable block gives a positive sign. The third block, corresponding to 
the neutral block, gives approximately 



exp 




1 



1 



-Ai^^e + e^ -AiAi^^-ak' 



{s)ds 



exp (-o-A;Vr^^) - 1 



Notice each of the three subsystems gives real value, since they clearly have real coefBcients. 
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and hence has the opposite sign as the dispersion parameter a. Notice, here we have used 
heavily the fact that the periodicity of the background solution implies 



and hence yields cancelation in averaging of the lower order terms. Combining these results, 
we find that the periodic Evans function has sign sgn(cr) for k ^ and /x ^ 1 as claimed. □ 

Notice that one could rework the above proof in the case where the background solution 
corresponds to a homoclinic orbit of the traveling wave ODE (12. ip . As a result. Lemma 13.21 
still holds in the solitary wave setting. This result seems to be new to the literature, and 
may give valuable insight to the transverse instability analysis in the solitary wave case. 

3.2 Low Frequency Analysis and Instability 

Now that we have a handle on the limiting sign of D{fj,,k,l) as /i — )• ±00, we turn our 
attention to determining the sign of the quantity D{0, k, 1) for |/c| ^ 1. For, once we have 
this information we see that the negativity of the orientation index 



provides a sufficient condition for transverse instability of the underlying periodic wave u. 
To this end, we utilize the methods of [J2j . which builds off the methods of |Br J] . to derive 
an asymptotic expansion of the function D{0,k,l) for \k\ <^ 1. However, it should be 
pointed out that the analysis in this case is seemingly more delicate than that considered 
by [J2] due to the fact that the unperturbed spectral problem at = 0, i.e. 



which does not define a Hamiltonian equation. In particular, this equation is not reducible 
to quadrature and hence one can not use integrability of the corresponding traveling wave 



{-Uxxx - f{u)x + CUx)x = 

to construct a basis for the monodromy operator at /i = 0, A; = 0. However, one can use 
the integrability of the traveling wave ODE (12. ip for the gKdV equation to construcllll three 
linearly independent solutions of equation (j3.5p : namely, the functions Ux, Ua, and ue are 
easily seen to satisfy 

C[u]ux = 0, C[u]ue = 0, C[u]ua = —1 
(see [BrJj for details) and hence provide a basis of solutions for the differential equation 

dxC[u]v = 0. 

^ Notice that here we are considering the traveling wave ODE as a formal differential equation with out 
any reference to boundary conditions. 




a-sgn{D{0,k,l)) 



(3.5) 



dl£[u]v = 0, 



ODE 
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Thus, we are missing one null direction of the (formal) operator d'^C[u]: we need a solution 
to the equation C[u]v = x. However, notice that the functions Ux and ue provide two 
Imearly independent solutions of the differential equatior^ C[u]v = 0, and hence one can 
use variation of parameters to solve the nonhomogeneous equation. After a straightforward 
calculation it is seen that 



sus{s)ds I ue{x) 



(x) :- 



is linearly independent from Ua^ 

and Ue and satisfies C[u](f) = x, and hence we can use 
these four functions to construct the corresponding monodromy matrix at the origin. Using 
perturbation theory then, we should be able to determine how these functions bifurcate as 
k varies but remains very small and hence be able to determine a leading order expansion 
of the periodic Evans function near k = 0. This is the content of the following lemma. 

Lemma 3.4. The following asymptotic relation holds in a neighborhood of k = 0: 
D{0, k,l) = - {PT - M2) {T, M}a,E [ak^y + 0{k^). 

Proof. We begin by writing the linearized equation ()3.2p as a first order system with coef- 
ficient matrix 



H(x, fi, k) 










V -ae - f"'{u)ul - f"{u)ux 



1 \ 

1 

1 

-2f"{u)ux - II -f'{u) + c J 



We now define the matrix W(x, fi, k) as the matrix solution of the first order system 
Y' = Yl{x, fi, k)Y such that 



(3.6) 



W(x,0,0) 



/ Ux 

u. 



■XX 

Uxxx 



Ua 
Uax 
Uaxx 
Uaxxx 



UE 
UEx 
UExx 
UExxx 



\ 



H^x 

^xx 
^xxx 



and we fix the initial condition W(0, ^, k) = W(0, 0, 0) for all (/i, /c) G C x M. Defining 
(5W(/x, A;) = W(x, )U, A;)|^^Q, a straightforward calculation gives 



5W(0,0) = 

\0 -V'{u^)V"{u-)Ta -V'{u-)V"{u-)Te V"{u^)V'{u^) xuE{x)dx J 



^ !q xux{x)dx \ 



/ 

V'{u_)Ta V'{u^)Te 

-T + V"{u-)^ xux{x)dx 

T 



-V'{u^) Jq xuE{x)dx 



Here, again, we consider the formal operator C[u\ without any reference to boundary conditions 
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Since 5 W(0, k) is analytic in k'^ by Lemma |3. 11 it follows that 

_ det(JW(0,fc)) _ 
^^"'^'^^-det(W(0,0,l)) " ^ ^ 

for <^ 1. Moreover, since the first column of the matrix (5W(0, /c) is 0{k'^), one can 
also easily see that the 0{k'^) variation in the ^ direction contributes to terms in D(0, fc, 1) 
of order 0{k^) near A; = 0. Thus, in order to compute the 0{k^) variation in the Evans 
function, we need only compute the 0{k'^) variation in the Ua, and ue direction^. 

Computing the necessary variations can be done using the variation of parameters for- 
mula. To this end, we define the vector solutions corresponding to Ux, Ua, ue and cp be 
given by Yi, Y2, Y3 and Y4, respectively and define 



—YjiT,o,k)\^^^=wiT,o,o) w(x,o,o)-i(y,- 



[x) ■ Cj) e4,dx, 





-V'{u^) 


V y"(n_)y'(n_) 



9m _ 
da 





~dE' 





where • represents the standard inner product on M^, ej is the j''^ column of the identity 
matrix on M^, and W(r, 0, 0) = W(0, 0, 0) + 5 W(0, 0) where 

/ 

W(0,0,0) = 

u L — V 



Moreover, in Appendix B we will show that 

/ - lo lo '^Eiz)dzds \ 

(3.7) W(x,0,0)~^e4 = 



\ 




-1 / 



V 



— X 

lo u{s)ds 
-1 



J 



which allows for straight forward computation of the necessary variations. Thus, we have 
5W(0,A;) = <5W(0,0)- | ^ Ay.(T, 0, ® e,- 1 ak^ + 0{k' 



and since det( W(0, 0, 0)) = 1 it follows that 



d 



D(0,A:,1) = det I 5W(0,0) - | ^ afc2^^(^' ^)L=o ® ) ) 



®Note that finding a useful expression for tfie C'(fc^) variation in tfie 4> direction would be a very daunting 
task, since one would have to apply variation of parameters to solve a problem with non-homogeneity which 
was constructed using variation of parameters. 
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Elementary row operations allow one to simplify the matrix involved in the above de- 
terminant. Indeed, a large simplification is possible if replaces the third row with the result 
from multiplying the first row by V"{u-) and adding it to the third row, and similarly re- 
place the fourth row with V"{u-) times the second row plus the fourth row. In particular, 
these operations imply the fourth entry in the first column is of order 0{k'^) and hence does 
not enter into the calculation. The resulting expressions are still too long to explicitly write 
out here, but can be easily handled by using a computer algebra system. Indeed, a direct 
calculation yields 

D{0, k, 1) =(^Q,^ + Q2^) V'{u^){T, M},,£;(afc2)2 + 

where Qi = (M^ - PT - MTu. + T^u^) and Q2 = {MT - T^u^). Finally, by noticing 
from ([22]) that 

l^'(n_)|^ = l and V'{u^)^ = u. 

we have 

*^+*^)''>-' = *^^-^^- 

which completes the proof. □ 

Remark 3.5. In WrJ^ . Bronski and Johnson derive the low-frequency expansion for the 
Evans function Dg^dvil^'i i"^) for the gKdV equation ()1.3p in the spectral parameter fi. In 
particular, it was shown that 

O..w(.,l) = 4d«(^(™),3^0(|,.|.) 

and hence an index which detects exponential instabilities to co-periodic perturbations was 
derived by comparing the sign of the above Jacobian determinant to the limiting behavior for 
/i — )• +00. We believe that it would be interesting and beneficial to derive the correspond- 
ing low-frequency expansion for the Evans function D{fi, k, k) considered in this paper, as 
it may give a better understanding of the connection (if there is any) between the genera- 
tion of transverse instabilities and uni- directional instabilities. Indeed, such analysis could 
serve as a starting point for a corresponding transverse stability analysis of periodic waves. 
While it seems natural that the above three-by-three Jacobian determinant should control the 
(spectral) low-frequency limit, it is not clear how the degeneracy of the gKP equation (cor- 
responding to the extra x-spatial derivative) affects the situation. In particular, the required 
calculations seem to be considerably more complicated than those considered in Lemma \3.4\ 
and a useful identification of the leading order behavior for <C 1 is yet to be obtained. 

Combining Lemmas 13.21 and l3. 41 we have a sufficient condition for exponential instability 
of a periodic traveling wave of the gKdV to long wavelength transverse perturbations in the 
gKP equation. This is the content of our main theorem, which we now state. 
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Theorem 1. Let u = u{-;a, E,c) be a periodic traveling wave of the gKdV equation (jl.3p . 
Then u is spectrally unstable to transverse perturbations in the gKP equation ()1.4p if the 
product a ■ {T, M}a.E is positive. 

Proof. Clearly the negativity of the orientation index 

sgn(Z)(0,A;, 1)) • lim sgn fc, 1)) 

/I— >- + oo 

for < ^ 1 implies the desired instability. Since PT — > by Jensen's inequality, 
the result follows by Lemmas 13.21 and 13. 4i □ 

Remark 3.6. Notice the instability detected in Theorem [7] is that of spectral transverse 
instabilities to perturbations which are co-periodic in the x-direction with low-frequency os- 
cillations in the transverse (y) direction. In the solitary wave case, a general criterion for 
spectral transverse instability was recently provided by Rousset ant Tzvetkov IRTl^ . While 
it seems plausible that such a criterion may exist when the underlying wave is spatially pe- 
riodic, we must note that the detected instability would not be to low -frequency oscillations 
in the transverse direction. Indeed, in JRTlf the transverse frequency must be large enough 
that the kernel of a particular linear operator be simple. Once such a frequency ^ is 
found, an implicit function type argument is used to prove the existence of a spectral curve 
H — )• {k{n),iJ,) in a neighborhood of fi = with (A;(0),0) = (A;o,0), thus proving spectral 
instability. In the present work however, we are proving the existence of a map k — )• fj,{k) 
defined for \k\ <^ 1 such that ^i{k) bounded away from zero for all k ^ and such that fj,{k) 
is an eigenvalue of (13. 2p for the associated small transverse frequency k. 

As previously noted Theorem [T] is only of interest when the underlying periodic wave is 
a stable solution of the corresponding one dimensional problem. In [BrJKj . the nonlinear 
(orbital) stability of such solutions to periodic perturbations in the gKdV equation was 
studied and the Jacobian {T, M}a,E was seen to play a significant role. Therein, it was 
shown that such waves can be nonlinearly stable to such perturbations regardless of the 
sign of {T, M}a,E, assuming certain conditions on the perturbation and other geometric 
quantities related to the underlying wave and the conserved quantities of the PDE flow. 
Theorem [T] however demonstrates the direct influence the sign of this Jacobian determinant 
has on the stability to perturbations in higher dimensional models (see also the recent 
work of Johnson [J2] concerning the transverse instability of periodic gKdV waves in the 
generalized Zakharov-Kuznetsov equations where the same Jacobian was seen to control the 
low-frequency behavior of the corresponding Evans function). In particular, we immediately 
have the following interesting (and seemingly unexpected) corollary. 

Corollary 1. A periodic traveling wave of the gKdV for which {T,M}a,E 7^ can never be 
spectrally stable to transverse perturbations in the gKP equation for both signs of dispersion. 

Remark 3.7. It should be noted that, unlike the ODE case, there are no general theorems 
insuring that spectral instability implies nonlinear instability. However, in the recent work 
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JRTSf it was shown in the solitary wave context that, indeed, spectral transverse instability 
of KdV waves in the KP-I equation (as described in the introduction) can be converted to 
a nonlinear instability result. It seems plausible that the methods utilized could apply in 
our case in order to convert Theorem{l\ into a nonlinear instability result (at least in the 
case of periodic transverse perturbations). This would be an interesting direction for future 
investigation. 

We now point out several corollaries of Theorem[TJ In the case of the KdV equation (II. ip . 
it is known that all periodic traveling wave solutions are both spectrally stable to localized 
perturbations (see [BDj ) and nonlinearly (orbitally) stable to co-periodic perturbations 
(see [BrJKj or |Jlj). In particular, the transverse stability of such solutions in the KP 
equation (11.2p is of interest in this case. Moreover, it was shown in |BrJK| that the Jacobian 
{T, M}a,E can be expressed as 



where disc(i?(-)) represents the discriminant of the polynomial R. Since V' is clearly strictly 
convex in this case, it follows by Jensen's inequality that 



For an alternate proof of this fact, see [J2j . Moreover, notice that for any {a,E,c) G M'^ 
for which (j2.ip admits a periodic solution of the KdV the equation E = V{u; a, c) has 
three solutions in u and hence the discriminant must be positive. Therefore, we have that 
{T, M}a,E > for all periodic traveling wave solutions of the KdV equation. This proves 
the following corollary of Theorem [H 

Corollary 2. All periodic traveling wave solutions of the KdV are unstable to long wave- 
length transverse perturbations in the KP equation when o" > 0. 

It is interesting to note that it is known that all solitary wave solutions of the KdV 
are transversely unstable in the KP equation when a > 0. Thus, Corollary [2] seems to be 
somewhat expected. Moreover, it turns out that the Galilean invariance of the KdV implies 
we can always choose a = 0, and hence (up to translation) the periodic traveling waves 
of the KdV form only a two parameter family of solutions. This family can be expressed 
explicitly in terms of the Jacobi elliptic function as 



where uq is an arbitrary parameter (taking the role of E) and k is the elliptic modulus. 
We refer to such a solution as a cnoidal wave solution of the KdV. As a result of Corollary 
[21 it follows that all cnoidal wave solutions of the KdV are unstable to long wavelength 
transverse perturbations in the KP-I equation. 





12 disc{E -V{-,a,c 




u{x, t) =uo + Uk'^K^ cn^ {k {x + {Sk'^K^ - Ak^ + uo) t) ,k) , 
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We now move on to consider a periodic traveling wave solutions of the focusing mKdV 
equation 

Ul — Uxxx ~t~ U tlx 

with positive wave speed c > 0. When a = 0, the corresponding traveling wave ODE admits 
two distinct classes of periodic solutions: when E > the wave can again be expressed in 
terms of the Jacobi elliptic function cn, while when E < there either exists no solution 
(if l-E"! is sufficiently large) or the solution can be expressed in terms of the Jacobi elliptic 
function dn and hence represents a dnoidal wave. In the recent work of [Br JKj . it was shown 
that both sets of solutions are nonlinearly (orbitally) stable to co-periodic perturbations, 
although the cnoidal solutions of sufficiently long wavelength were shown to be (spectrally) 
unstable to periodic perturbations of large periocll^ (see also |DK| ) . Moreover, the sign of 
the Jacobian {T, M|„ k was analyzed for both the cnoidal and dnoidal wave solutions of 
the focusing mKdM^H and was seen to be positive for all dnoidal solutions and negative for 
all cnoidal solutions. As a result, we have the following corollary of Theorem [1] 

Corollary 3. AH cnoidal wave solutions of the focusing mKdV equation are unstable to 
long wavelength transverse perturbations in the focusing mKP equation with o" < 0, while 
the dnoidal wave solutions are unstable to such perturbations when o" > 0. 

Continuing, one can use the general elliptic function calculations of |BrJK| in the case 
of a power-law nonlinearity f{u) = u^^^ , p G N, to determine the sign of the Jacobian 
{r, M}a,E for any periodic traveling wave solution of (jl.3p in terms of moments of the 
background solution u with respect to the density 

1 

^E-V{-^' 

As such. Theorem [1] can be utilized to provide a transverse instability result for such power- 
law nonlinearities. In other cases, it seems that one must (in general) resort to numerical 
methods to approximate {T,M}a,E- 

4 Conclusions &; Discussion 

In this paper, we analyzed the spectral instability of a periodic traveling wave solutions of 
the generalized Korteweg- de Vries equation to long wavelength transverse perturbations in 
the generalized Kadomtsev-Petviashvili equation. In particular, we constructed a seemingly 
nonstandard orientation index by comparing the low and high frequency behavior of the 
periodic Evans function when the transverse wave number k is non-zero. We found that 

^"However, cnoidal waves of smaller period seem spectrally stable to perturbations of sufficiently large 
period. 

^^In fact, the situation was analyzed without the restriction of a = in which case all periodic traveling 
wave solutions of the focusing mKdV can not be expressed simply in terms of a Jacobi elliptic function. 
However, we only consider the case a = Q here for simplicity. 
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in the high frequency hmit, the Evans function -D(/U, k, 1) converged to zero as — )■ ±00, 
which is insufficient to conclude a instabihty theory. However, after taking into account 
higher order effects, it was found that by the periodicity of the underlying wave and the 
resulting cancelation in averaging procedures that the Evans function for non-zero transverse 
wave numbers favors a particular sign as — )• ±00 which is determined precisely by the 
dispersion parameter a: such a phenomenon seems to be new in the literature. Thus, an 
instability index follows by comparing the sign of a with the value D{0,k,l) for k ^ 0. 
Utilizing the methods of |BrJ ] and (J2J then, we were able to explicitly compute the leading 
order variation of the function D{0,k,l) in /c in terms of a Jacobian from the traveling wave 
parameters to the period and mass of the background solution. This Jacobian was shown to 
be (generically) non-zero in the physically important cases of the KdV and mKdV equations, 
and their resulting signs were inferred from the recent work of [BrJK] immediately yielding 
instability results in these cases. 

It is interesting to note that the Jacobian arising in the low frequency expansion of 
the periodic Evans function has already been seen to hold vital information concerning the 
stability of the periodic traveling wave solutions of the gKdV. Indeed, in [BrJKj and |J1] 
this Jacobian arose naturally in the nonlinear stability analysis of such solutions to periodic 
perturbations in the gKdV equation, while in |J2] it was seen again to control the low 
frequency behavior of the periodic Evans function when considering the spectral instability 
of a periodic gKdV wave to long wavelength transverse perturbations in the generalized 
Zakharov-Kuznetsov equation (which also arises in plasma physics). Thus, it would be very 
interesting to better understand the physical mGcining of tliG Jcicobicin -^^a^E 

as this 

may better illuminate the stability theories described above. 

Another interesting direction would be to complement the transverse instability analysis 
in this paper with a corresponding stability theory. That is, to derive sufficient conditions 
to guarantee a periodic traveling wave solution of the gKdV is transversely stable to per- 
turbations in the gKP equation. While this certainly may be possible using the Evans 
function techniques of this paper, much more delicate analysis is needed. In the context of 
the KdV equation, one may be able to use the integrable structure of the KP-II equation 
to prove transverse spectral stability. Such techniques are prevalent in the solitary wave 
theory and have recently been employed in the periodic wave setting to prove spectral sta- 
bility to localized perturbations for several model equations (see [BDJ . [BDNJ, and [NBJ). 
When this integrable structure does not exist, however, it may be more natural to consider 
a variational characterization of the stability problem and extend the methods of |GSS] . 
[JT], and pTjiq. 

We also note that, after the completion of this work, our attention was brought to the 
recent work of Rousset and Tzvetkov [RT2j in which the authors considered the transverse 
spectral instability of KdV solitary waves to perturbations in the KP-I equation. In contrast 
to the ODE techniques utilized in this paper, the authors use variational techniques to 
present a rather elegant and simple approach relying only on properties of the differential 
operators involved which are rather easy to check (due to the self-adjointness of the operators 
involved). We believe it would be interesting as a future direction of study to see if these 
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techniques could apply in the case where the underlying wave is spatially periodic and to 
compare the results to those derived in this paper. Such a comparison would hopefully help 
illuminate the mechanism behind the instability. On the other hand, the ODE techniques 
used in this paper are quite robust allowing for straightforward numerical implementation 
(plotting Evans curves and computing winding numbers) and applying not only to equations 
with a Hamiltonian like structure but also to more complicated situations arising in the 
context of systems of nonlinear conservation laws. Moreover, the degeneracy in the high 
frequency analysis is of independent interest adding new techniques to the tool box in the 
study of nonlinear dispersive waves utilizing asymptotic tracking/reduction results familiar 
from shock wave analysis in the conservation law setting. 

Finally, as pointed out in the text the high frequency analysis conducted in this paper 
translates directly to the solitary wave setting, and hence opens the door to an analogous 
instability theory using Evans function techniques. In fact, the high frequency limit seems 
easier to discern in the solitary wave case as one does not have to deal with averaging 
effects of the coefficient functions. While we have not yet conducted the relevant low 
frequency analysis, this could provide valuable insights concerning the transverse instability 
of a solitary traveling wave of the gKdV in the gKP equation. 

Acknowledgement. Thanks to Jared Bronski, Bjorn Sandstede, Todd Kapitula, and 
Bernard Deconinck for many useful conversations regarding the early stages of this work. 
Their comments and suggestions were invaluable to this project. We would also like to 
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A A Block- Triangular Tracking Lemma 

Consider an approximately block-triangular system 

(A.l) W' = AP{x)W := ^f^ M^^^P)^^ 

where O is a uniformly bounded matrix, 5{x) scalar, and p a vector of parameters, satisfying 
a pointwise spectral gap condition 

(A.2) mincj(5RMi) - maxo-(KM2) > ??(x) > Ofor all x. 

(Here as usual ^N := (1/2)(A^-|-A^*) denotes the "real", or symmetric part of A^.) Then, we 
have the following block-triangular version of the tracking/reduction lemma of |MaZ31 IPZ] . 
For related results, see [HLZj . 

Lemma A.l. Consider a system (jA.ip under the gap assumption ()A.2p . with Q uniformly 
bounded andrj G L^^^. If snp{6/r]){x) is sufficiently small, then there exists a unique bounded 
linear transformation 

(A^3) S^(^»), 
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possessing the same regularity with respect to p as do coefficients Mj and N , such that the 
change of coordinates W = SZ converts the approximately triangular system ()A.ip to an 
exactly block triangular system 

(A.4) Z' = A^{x)Z := {^^^ ^y^'^P^^ 

where 

(A.5) Ml := Ml + $iV, M2 := M2 - ^N, N := iV, 

with 

sup < C sup(5/77) 

and 

(A.6) \i^{x)\<C e^>^'^'^'5{y)dy, \^{x)\ < C f 6{y)dy. 

J X J —00 

where the constant C depends only on the size of Q . 

Proof. By the change of coordinates x — >• x, 5 — )■ 5 := 5/r} with dx/dx = ri{x), we may 
reduce to the case rj = constant = 1 treated in |MaZ3j . Dropping tildes, we find by direct 
computation that Z = S~^W satisfies ()A.4p for S of form ()A.3P if and only if (jA.SP and 

^>' = (M2$-$Mi) + Q(^>), 

where Q is the quadratic matrix polynomial Q{^) ■= 5@ — $A^$. Viewed as a vector 
equation, this has the form 

with linear operator := M2<5 — 'I'Mi. Note that a basis of solutions of the decoupled 
equation = A^<I> may be obtained as the tensor product ^ = (pep* of bases of solutions of 
(j)' = M2(j) and cp' = —M^cp, whence we obtain from (jA.2p 

(A.7) e^"<C7e-''^ for z > 0, 

or uniform exponentially decay in the forward direction. 

Thus, assuming only that <I> is bounded at —00, we obtain by Duhamel's principle the 
integral fixed-point equation 

(A.8) $(x) = r$(x):= r e^(^-f)Q(^>)(y)dy. 

J —00 

Using (IA.7p . we find that T is a contraction of order 0{5/r}) for <^ on a ball of radius of 
the same order, hence (jA.SP determines a unique solution for 5/r} sufficiently small, which, 
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moreover, is order 5/rj as claimed. (Here, we are using the normalization rj = 1.) Finally, 
substituting Q($) = 0{6 + |^>p) = 0{6) in (IXSll . we obtain 

|$(x)| <C r e''^''-y'^ 6{y) dy 

J —oo 

in X coordinates, or, in the original x-coordinates, ()A.6p . Regularity with respect to param- 
eters is inherited as usual through the fixed-point construction via the Implicit Function 
Theorem. □ 

Remark A.l. Although we do not use it here, an important observation of the above proof is 
that for rj constant and 8 decaying at exponential rate strictly slower that e"^^' as x ^ +oo, 
we find from ()A.6P that <I>(x) decays like 5/ri as x ^ +oo, while if 6{x) merely decays 
monotonically as x ^ — oo, we find that ^(x) decays like {5/r]) as x ^ —oo. 

Remark A. 2. Though we do not use it here, an important observation of JMaZ3[ \PZf is 
that hypothesis (IA.2[) of Lemma \A.1\ may be weakened to 

min(T(3f?Mf) - max a($ftM|) > r/(x) + a{x,p) > 

with no change in the conclusions, for any a. satisfying a uniform L^ bound \a{-,p)\ii < Ci. 
(Substitute e'^^ < Ce^^e~^^ for (lA.Tp . with no other change in the proof.) This allow us to 
neglect commutator terms in some of the more delicate applications of tracking: for example, 
the high-frequency analysis of fMaZS^ . 

Remark A. 3. In the special case that A^ is T-periodic in x, we obtain by uniqueness that 
<I> is T-periodic in x as well. 



B Variation of Parameters Calculation 

In this appendix our goal is to justify equation p.7p . which was seen to be a crucial step 
in the low- frequency analysis of Section 3. For brevity, however, we will consider only the 
most difficult case. To begin, define the scalar valued function ^0(2;) := e\ W(x, 0, 0)~^e4, 
where W(x, 0, 0) is defined in (j3.6p and f represents the vector adjoint, and note from ()3.7p 
we wish to prove Aq = — UE{z)dzds. By definition, we see that 

^0 = {UaxxUEx — UaxUExx) <t> + {Ua^Exx — UaxxUE) 4>x + {u 

ax'^E '^afJ'Ex) (Pxx- 

Using (12. 2p . we have that 

UaUExx - UaxxU - E = -UaV" {u)ue + Ue {V" {u)Ua - l) 
= -UE- 

Noting that {UaUEx - UaxUE) = UaUExx " UaxxUE, it folloWS that 

rx 

UaxUE - UaUEx = / UE{s)ds. 

Jo 
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Moreover, using (j2.2p along with the fact that 

Ux {UaxUE - UaUEx) = UUe - Ua, 

it follows that 

Ux {UaxUExx " UaxxUEx) = -UExUx - V"{u) {uUe - Ua) 

= -UExUx - V"{u)Ux {UaxUE " UaUEx) 



= -UExUx - V"{u)Uj; / UEis)ds. 

Jo 

Therefore, using the above equalities along with the definition of the function (j), we have 

J-X 

Aq = —X / UE{s)ds 
Jo 

+ (J^Ex + V"{u) j UE{s)ds^ Ux - UEUxx + Uxxx j UE{s)ds^ j SUE{s)ds 

+ (^UEx + y"{u) j UE{s)ds^ UE - UEUEx + UExx j UE{s)ds^ j SUs{s)ds. 



Now, as above one can show 



/•X \ rx 

UEx + V"{u) J UE{s)ds\ Ux - UEUxx + Uxxx j UE{s)ds 

px 

= UExUx - UEUxx + {V"{u)Ux + Uxxx) / UE{s)ds 

Jo 



= 1 

—V"{u)ux by Similarly, it follows that 



and hence 



rx \ j-x 

UEx + V"{u) j UE{s)dsj UE - UEUEx + UExx j UE{s)ds = 

PX PX 

Aq{x) = —X / UE{s)ds + I suE{s)ds 
Jo Jo 

(s — x)uE{s)ds 
UE{z)dzds 



JO 



as claimed. The rest of the derivation of equation (j3.7p is handled similarly, although the 
necessary calculations are considerably simpler. 
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